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ABSTRACT 



Motivated by the significant overdensity of background bright quasars 
recently detected behind the foreground clusters of galaxies on scale of 10 
arcminutes, we have investigated the possibility of attributing the quasar-cluster 
associations to gravitational lensing by large-scale matter inhomogeneities. 
Based on the conventional lensing models, we have shown that the reported 
quasar overdensity is unlikely to be generated by cluster matter alone. The 
situation does not change even if all the clusters of galaxies which follow their 
spatial two-point correlation function arc taken into account, while matter 
clumps on scale of > 20 Mpc are also found to be unable to provide the required 
mass surface density since their density contrast is strictly limited by the 
anisotropy measurements of the cosmic background radiation. Moreover, we 
have pointed out that the influence of a nonzero cosmological constant on the 
quasar-cluster associations is very minor. We conclude that either the observed 
quasar number counts have been seriously contaminated by the magnification 
bias of matter inhomogeneities of the universe or there should exists some 
intercluster matter on scale of less than ~ 20 Mpc, e.g. from cluster-galaxy 
correlation, whose mean cosmic density is about an order of magnitude higher 
than that of clusters of galaxies. 



Subject headings: cosmology: gravitational lensing — galaxies: clusters: general 
— large-scale structure of universe 
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1. Introduction 

A significant overdensity of background bright optical/radio quasars have been recently 
detected on scale of 10' around foreground Zwicky/Abell clusters (Rodrigues-WiUiams & 
Hogan 1994; Wu & Han 1995; Rodrigues-WiUiams & Hawkins 1995; Seitz & Schneider 
1995). Although these unusual quasar-cluster associations are generally believed to be the 
result of statistical lensing of quasars by foreground gravitational potential, cluster matter 
alone is far from explaining the observed amplitudes of the quasar overdensity behind 
clusters. It then seems that the matter inhomogeneities on even larger scale (> 10 Mpc) 
traced by galaxy clusters should be taken into account in the explanation of the reported 
quasar-cluster associations. 

While the overdensity of background quasars on the similar scale around foreground 
galaxies was detected a few year ago (Fugmann 1988; 1990; Bartelmann & Schneider 
1993b;1994), the large-scale matter clumps that galaxies are associated were also advocated 
in order to produce the quasar-galaxy associations. Using N-body simulations of galaxies 
formation, Bartelmann & Schneider (1993a) did find a correlation of high redshift quasars 
with low redshift galaxies in the scenario of magnification bias by the matter of galaxies 
and their surrounding large scale structures. Their results indicate that galaxies, and 
probably clusters of galaxies, contribute a minor effect on the quasar overdensity on scale 
of arcminutes. Interestingly, the recent work by Wu, Zhu & Fang (1995) shows that even 
on small scale of arcseconds the quasar-galaxy associations are actually generated mainly 
by the cluster matter rather than the galaxies. 

Compared to the quasar-galaxy associations, the quasar-cluster associations deal with 
the matter distributions on scale of ranging from ~ 1 to ~ 10 Mpc, on which galaxy 
clusters are strongly correlated, revealed by their two-point correlation function ^(r). It 
is timely and necessary to address the following question: Is galaxy cluster clustering 
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described by ^ (r) able to provide enough gravitational matter to act as lens for the reported 
quasar-cluster associations ? A definite answer to such a question today relies on the 
numerical study of various models of formation of large-scale structure of the universe. 
However, an analytic investigation, as we will make in this letter, of the quasar-cluster 
associations in the scenario of gravitational lensing by various matter clumps may supply 
us with a very useful clue to the matter distribution on large-scale of the universe. 



2. Quasar enhancement factor 

The overdensity of background quasars at an angular distant 9 around a foreground 
galaxy cluster is described by the enhancement factor q{9) (Narayan 1989) 

iV[<m + 2.51ogM^)] 1 _ N[>S/m] 1 

' N{<m) fi{e) N{>s) ix{ey ^ ' 

where N are the intrinsic quasar number counts above a limiting magnitude m or a fiux 
threshold S and ^{9) is the lensing magnification introduced by the foreground gravitational 
potential. This equation accounts for both the magnification effect (2.5 log ;U or S/ ^) and 
the area distortion due to light deflection by foreground matter. To compare with 

the measurements of quasar-cluster associations that search for quasar number excess over 
a range with radius 9 around the cluster center, the average enhancement factor {q{9)) is 
employed: {q{9)) = 2 q{9)9d9/9^. 

If we assume that the observed quasar number counts are not significantly 
contaminated by gravitational lensing due to the matter clumps in the universe, then 
the number-magnitude relation A^(< B) from Boyle, Shanks & Peterson (1988) and the 
source counts A^(> S) at 5 GHz from Langston et al. (1990) can be adopted for the 
optically-selected and the radio-selected quasars, respectively. However, it should be noted 
that the radio counts N{> S) contain both quasars and galaxies and the fraction of quasars 
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in A^(> 5") varies with the flux threshold. Therefore, the employment of A'"(> S) in the 
study of quasar-cluster associations can only be regarded as an approximate estimate of (q) . 
Furthermore, quasars are treated as pointlike sources, which would be suitable for clusters 
and large-scale matter clumps as lenses. 

The observational data of the four searches for quasar- cluster associations provide 
actually the variations of (q) with the search distance and/or the limiting magnitude. 
We don't intend to flt the curves from the theoretical modeling of lensing systems as did 
by Wu & Han (1995) and Rodrigues- Williams & Hawkins (1995). Instead, wc adopt the 
one significant result of (g) measured at a fixed 9 and a limiting magnitude (or fiux). 
Rodrigues- WiUiams & Hogan (1994) have exphcitly given the enhancement factor within 
6 Zwicky radii and B < 18.5. The enhancement of (g) versus 9 becomes nearly unity at 
6' > 5' in Wu & Han (1995), and therefore, 9 = 4' seems to be a reasonable "edge" for 
their sample. Rodrigues- WiUiams & Hawkins (1995) show the variations of (q) against the 
quasar hmiting magnitude for 9 ^ 0.12". We take the most significant value at S = 18. 
Scitz & Schneider (1995) choose the quasars from the 1 Jy radio source sample at 5 GHz 
but actually use the optically identified sources which are mostly quasars or BL Lac objects 
at high redshift Zg > 0.5. Their most significant signal appears at S < 19 and 1. We 
adopt the value of (g) fa 1.3 at 9 Zwicky radii. In particular, we keep both the hmiting 
magnitude and the radio fiux threshold, which helps to finger out the reliability of using 
radio source counts A'"(> S) in the evaluation of (q). Unfortunately, we cannot read out the 
uncertainties from their data but learn that the result has a very high significance of up to 
98%. Table 1 summarizes these four measurements, in which the errorbars in {q)obs are the 
laVN errors arising from the estimates of both quasar density over the association area 
and the mean quasar density. 
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3. Modeling of quasar overdensity behind clusters 

Now we work with the lensing models of the matter inhomogeneities associated with the 
foreground clusters and test what would be required to explain the reported quasar-cluster 
associations. We adopt a flat cosmological model with a cosmological constant Aq, i.e., 
f^o + Ao = 1, and Hq — 50 h^o km s~^ Mpc~^. 

Clusters of galaxies were naturally thought to be the deflectors for the observed 
quasar-cluster associations. We can evaluate the contribution to (q) from clusters by 
modeling the cluster matter as a singular isothermal sphere, which is characterized uniquely 
by its one-dimension velocity dispersion Cy. This profile is the simplest but more or less 
reasonable model for dark matter distribution in clusters of galaxies. We count both the 
primary quasar images and secondary ones, if any, which are gravitationally magnified by 
a factor of fj,{9) = |1 — 9E/d\^^, where the Einstein radius is 9e = 4:7r{av/c)'^Dds/ Ds, and 
we use Dd, Ds and D^s to denote the angular diameter distances to foreground lenses, to 
background sources and from lenses to sources, respectively. In Table 1 we hst the cluster 
velocity dispersion which is required to produce the observed enhancement for each of the 
four measurements. Apparently, the resulting cr^ is substantially larger than any realistic 
values for clusters of galaxies. Taking the mean cluster velocity dispersion as 1000 km s~^, 
we can estimate that the gravitational mass (M ~ cr^) responsible for the quasar-cluster 
associations is an order of magnitude higher than the presently known total cluster mass. 
Meanwhile, it is seen that the real matter clumps which generate the quasar- cluster 
associations must deviate from the distribution since one cannot use a single velocity 
dispersion parameter to reproduce all the observed {q) . 

If the large-scale matter inhomogeneities traced by clusters of galaxies contribute 
an additional mean surface mass density E to cluster matter, the Einstein radii 6e 
of the background quasars will be increased by a factor of (1 — S/Scrit)~^, where 
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Scrit = {c^ I^T^G){.Ds/ D^D^s) is the critical surface mass density (Turner, Ostriker & 
Gott 1984). The image magnification now reads ^{9) = |1 — 9E/d\^^{l — S/Scrit)^^- 
Therefore, if the large-scale matter inhomogeneities have a mass density comparable to the 
critical one, the magnification factor can be greatly enhanced. Alternatively, T^crit would 
be smaller in a cosmological constant dominated universe than in the matter-dominated 
universe, i.e., the same uniform matter sheet would act as more efficient lens in a Aq 
dominated universe (Wu et al. 1995). Table 1 gives the mean surface mass density S of 
the large-scale inhomogeneities that are needed to explain the quasar overdensity around 
clusters. Recall that the surface mass density at the cluster center with core radius of Tc is 
0.087(a^/103 km s"^)^ (r,/0.25 Mpc)-^ h^^ g cm and the minimum critical density for a 
source at = 2 is T^cnt = (0.41,0.28) h^o g cm~^ for VLq = (1.0,0.2). Thus, the large-scale 
matter clumps should have their surface mass density comparable to the one at the cluster 
center in order to act as the lenses for the observed quasar enhancement around clusters It 
is noticed that S deduced from the radio bright quasar associations with clusters is a factor 
of ~ 2 larger than the one from the optical quasar samples. As we have mentioned before, 
this is due to the contamination of radio galaxies in the radio source counts A^(> S) we 
adopted. The data of Seitz & Schneider (1995) illustrate very well this effect: The optical 
quasar number-magnitude relation A^(< B) results in a E that is a factor of about 2 smaller 
than the value given by the radio source counts A^(> S) for the same set of quasar data. 

EDITOR: PLACE TABLE | HERE. 

In summary, the mass distributions for clusters of galaxies fail in reproduction of 
the observed quasar-cluster associations, no matter how massive they would be. A uniform 
surface matter sheet is found to be an accepted model as long as its surface density reaches 
a value comparable to the one at the cluster center. The introduction of the cosmological 
constant does not significantly reduce the demand for such a high density inhomogeneity on 
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large-scale. 

4. Contributions of large-scale matter clumps 

We now discuss the mass contributions from the cluster-cluster correlation. The matter 
clustering that clusters of galaxies trace on large scale can be quantitatively described by 
the cluster spatial two-point correlation function ^(r) = (r/rcc)~^'^, where the correlation 
length is ~ 40 Mpc h^^ (Postman, Huchra & Geller 1992). Since ^(r) diverges at r = 0, 
we truncate ^(r) when r < tq. The probability of finding a cluster in the surface element 
2Tr(^d( at distance ( from a cluster on the plane perpendicular to the line of sight is 

dP{C) = AnnCdC + e(r)(l + z,Y]-j0L=, (2) 

in which n is the mean cluster number density and e accounts for the evolution of ^{r). The 
expected mass contribution from all the clusters following ^(r) can be computed by the 
integration of m{()dP{() over ( from to oo. Here m{() is the cluster mass within d( of (. 

We again adopt a singular isothermal sphere model for individual cluster with velocity 
dispersion a^. Moreover, the cluster matter distribution is truncated at the cluster 
gravitational radius Rc so that the cluster mass is Mc = 2a'^Rc/G. We consider only those 
excess population relative to the "background" cluster of mean density n, i.e., we take out 
the factor "1" in [1 -|- ^(r)]. The expected mean surface mass density over the area 
provided by clusters with mass Mc is [/ m{()dP{()] /{7r(^), which reads 

E(Co) = 4nMcrcc(l + zaY F(Co, Tcc, Rc, ro), (3) 

where 
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and 



^(ro, C) 



roo dx 1^ <r r • 

Jrn/C ^0.8, /^2_i ) S ^ '0) 



1.84, C > r-o. 

Summing up the contributions from all kinds of clusters with different mass and number 

density gives rise to the expected mass surface density: ^{Co) — 40cPo(l + ZdYrcc F, in 
which represents the fraction of the total cluster matter in the matter (with the critical 
mass density po) of the universe. For the typical cluster radii of i?c = 3 - 5 h'^Q Mpc and 
the smallest cluster separation of Tq = 5 - 10 h'^^ Mpc, numerical computation shows that 
F f=i 2 ~ 3 over the range of C = 1 - 20 h'^Q Mpc which is comparable to the search distances 
from cluster centers in the measurements of quasar-cluster associations. Therefore, 

E = 0.0117, f^^V' (^] /l5ogcm-^ (6) 



1.15 J V3. 

where we have assumed a stable clustering model e = —1.2 and converted the comoving 
surface mass density into a physical one by multiplying a factor of (1 + z^)^. These two 
factors do not significantly alter our following result since clusters involved are at relatively 
low redshift. It appears that even if we take flc = = ^, the cluster matter provided by the 
cluster-cluster correlation is of an order of magnitude lower than the surface mass density 
required to produce the observed overdensity of background quasars around foreground 
clusters. 

The mass surface density from matter inhomogeneities on scale of larger than the 
coherence length (~ 50 Mpc) of the cluster-cluster correlation can be estimated though 

fOO / R \ 

L ^''^'"^ " '""^ ~ '"^^^ " ^-^^ ( lOOMi^ ) ^ 

where 5 is the mean present density contrast over scale of R. However, the evaluation 
of 5 on scale of larger than ~ 10 h^^ Mpc is sharply constrained by the measurements 
of temperature fluctuation l\T jT of the cosmic background radiation on various angular 
scales. Using the simple model for a spherical density perturbation (Fang & Wu 1993), we 
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can set an upper limit on E in terms of the recent results of AT/T. It turns out that the 
resulting E from any mass clumps on scale of greater than ~ 20 h'^Q Mpc is at least an order 
of magnitude smaller than the mass surface density required to explain the quasar-cluster 
associations. Therefore, it is unlikely that the observed quasar-cluster associations can be 
attributed to the lensing effect by large-scale (> 20 h'^^ Mpc) structures of the universe. 



5. Discussion and conclusions 

We have shown that the strong associations of background quasars with foreground 
clusters on scale of ~ 10 arcminutes cannot be interpreted as the statistical lensing by 
clusters of galaxies. The situation does not improve even when all the cluster matter 
that follow the two-point cluster-cluster correlation function is involved, while the matter 
contribution from large-scale structures (> 20 h'^Q Mpc) is strongly constrained by the 
measurements of the temperature anisotropics of the cosmic background radiation. We are 
limited to very few possibilities to solve the puzzle of quasar-cluster associations. 

An intuitive speculation is that the reported quasar-cluster associations are statistical 
variations arising from the quasar/cluster selections. Rodrigues- William & Hogan (1994) 
and Seitz & Schneider (1995) have already pointed out that the patchy dust obscuration 
cannot explain their observations. Alternatively, the background quasar clustering is 
detected only at r < 60 h^^ Mpc (Mo & Fang 1993). This clustering scale is comparable 
with the angular separation in the quasar-cluster associations but is much smaller than 
the spatial separation of the selected quasars. Cluster-cluster autocorrelation seems to be 
another possibility. However, if the background quasars are detected randomly on the sky, 
there would be no angular correlation between quasars and clusters even if the clusters are 
auto-correlated, which has been shown by Rodrigues- WiUiams & Hogan (1994) using their 
data and also by our simulations. 
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Another way out of the difficulty is that there exists a large amount of unseen matter 
between clusters of galaxies on scale of less than ~ 20 Mpc. Recall that we did not include 
the contribution of the intercluster matter in the above discussion. It is hard to figure out 
the distribution of this dark matter, but it should be massive enough to provide a surface 
density of as high as that required by the quasar-cluster associations. We will employ the 
N-body simulation of formation of clusters and large-scale structures to further study the 
issue (Wu, Fang & Jing, 1996) 

It may be possible that the observed background quasar counts deviate significantly 
from their intrinsic ones. The fact that quasars are strongly associated with the foreground 
galaxies and clusters indicates that the observations may preferentially select those quasars 
whose angular positions appear to be close to the foreground matter clumps. Unfortunately, 
previous studies (Schneider 1987;1992; Pei 1995; references therein) about the magnification 
bias on the observed quasar number counts reached a controversial result, depending mainly 
on our current knowledge of the distribution of lensing objects in the universe. It deserves 
to be investigated whether or not the quasar counts have been seriously contaminated by 
the lensing effect due to large-scale matter inhomogeneities. Meanwhile, our theoretical 
prediction of q depends sensitively on the adopted quasar counts (Boyle et al. 1988), which 
may have large uncertainties. Recall that a different quasar number-magnitude relation is 
derived by Hawkins & Veron (1993). 

Finally, we have tested the possibihty of attributing the quasar-cluster associations to 
the cluster environmental effect from the gravitational matter of cluster-galaxy correlation 
and will present the result elsewhere (Wu et al. 1996b). 

We thank an anonymous referee for her/his valuable suggestions. WXP was supported 
by the National Science Foundation of China and a World-Laboratory fellowship. 
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Table 1. Quasar-cluster associations: observations and models. 



clusters quasars {zdY {zs)^ 6" {q)obs {(Ty/IO^Y S'^ ref 



Zwicky 


B < 18.5 


0.2 


1.8 


52 


1 7+0.5 


5.31};^ 




0.081°;°^ 


1 


Abell 


5'> 2 Jy 


0.1 


2.0 


24 


1 7+0.5 
^■'-0.5 




0.281°;!° 


0.251°;°^ 


2 


UKJ287S 


B < 18.5 


0.15 


1.5 


7.2 


'^■^-0.2 


9 O+0.2 


l2+° °2 

'-'•-'-'^-0.02 


11+001 


3 


Zwicky 


< 19 


0.2 


1 


78 


~ 1.3 


4.3 


0.06 


0.05 


4 



> 1 Jy 5.6 0.11 0.10 



REFERENCES. - (l)Rodrigues-Willianis & Hogan 1994; (2) Wu & Han 1995; (3) Rodrigues & 
Hawkins 1995; (4) Seitz & Schneider 1995. 
'^Mean cluster redshift 

''Mean quasar redshift 

'^Search range in arcminutes 

'^Required cluster velocity dispersion in units of 1000 km s~^ 
''Required surface mass density in g cm~^ h^o for = 1 and Aq = 
^Required surface mass density in g cm~^ h^o for flo = 0.2 and Aq = 0.8 
^Clusters in UKJ287 field 
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